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ABSTRACT 

Context. The first substellar companion of possibly planetary mass around a normal star, GQ Lup, has been directly imaged 
(Neuhauser et al., 2005). Besides the unknown formation history, the mass of such an object is a criterion to decide about its true 
nature. 

Aims. We aim to determine the physical properties of the GQ Lup companion - effective temperate (r cff ) and surface gravity (log g), 
and thus its mass independently from evolutionary models. 

Methods. We use the adaptive optics integral-field spectrograph SINFONI at the VLT for near-infrared spectroscopy from 1.1 to 2.5 
fjm with a resolution of R = 2500-4000. We compare these spectra with synthetic atmospheric models (GAIA v2.0 cond). 
Results. From the complete set of spectra we find a consistent effective temperature and surface gravity of T c ff = 2650 ± 100 K and 
log g = 3.7 ± 0.5 dex. Combined with a slightly revised luminosity of log L/Lq = -2.25 ± 0.24 dex for the companion, we determine 
a radius of R = Rj„ p and thus a mass of ~ 25M Jup . The uncertainty of this value is rather high. Due to the large uncertainty of 

the surface gravity, the mass could range from 4 to 155 Mj„ p . By comparing the paramaters of the companion of GQ Lup to the ones 
of 2MASS J05352 184-0546085, published by Stassun et al. (2006), we conclude that the companion to GQ Lup A has a mass lower 
than 36 ± 3 M Jup . 
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1. Introduction 



\ iNeuhauser et al.l (I2005I) and lMugrauer & NeuhauseJ d2005l) pre- 
sented a faint common proper motion companion to the 1 + 1 
Myr old T-Tauri star GQ Lup. Based on both, a low resolution 
NACO spectr um (R~700) in comparis on to early GAIA model 
atm ospheres dBrott & Hauschildtl 120051) and evolutionary mod- 
els (IWuchterlL 120051) . the companion was found to have spectral 
type M9 to L4 and a mass of up to a few Jupiter masses. 

The mass estimate was, however, model-dependent. Due to 
the young age of GQ Lup, the most widely us e d evo lutionary 
models from t he Tucson group dBurrows et all 1 1993b and the 
Lyon group dChabrier et al.L l20Q0h for l ow mass objects are not 
applicable, see, e.g., Baraffe et al. (2002|). 

The direct determination of the mass from effective temper- 
ature, surface gravity and luminosity was hampered by the low 
signal-to-noise ratio and low resolution of the NACO K-band 
spectrum that yielded T e ff and logg with wide error margins 
only. Moreover, the fit of the observed spectrum against a grid 
of synthetic spectra from GAIA dusty atmosphere models did 
not deliver consistent results for the CO band features an d the 
overall slope of the spectrum (see e.g. lGuenther et all (120051) ). 

It was thus desirable to obtain new near-infrared spectra with 
higher signal-to-noise ratio and higher resolution over a wider 
spectral range. Because of the large magnitude difference be- 
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tween GQ Lup A and its companion and their small separation of 
~0.7"the use of adaptive optics systems was mandatory. Hence, 
we used the adaptive optics integral-field spectrograph SINFONI 
at the Very Large Telescope (VLT) that provides a spectral res- 
olution of R-2500 in the J-band and R-4000 in the H and ri- 
band. 

We present the observations and describe the data reduction 
in section 2. The spectral synthesis is outlined in section 3. We 
present the results of the fitting process in section 4, where we 
end with conclusions and a short discussion about the nature of 
the GQ Lup companion. 

2. Observations and data reduction 

2.1. Observations 

SINFONI (Spectrograph for INfrared Field Observations) is 
a combination of two instruments. First, a MACAO (Multi- 
Application Curvature Adaptive Optics) type curvature adaptive 
optics (AO) module with visual wavefront sensor and 60 actu- 
ators. Second, a mid-resolution near-infrared spectrograph with 
an integral field uni t: SPIFFI (Spectrograph for I nfrared Faint 
Field Imaging), see lEise nhauer et al. (2003J) and iBonnet et al.l 
(2004). The advantage of this instrument that combines the spa- 
tial resolving power of an 8m telescope with a mid-resolution 
spectrograph is the absence of wavelength dependent slit losses 
that occur on normal spectrographs with a narrow entrance slit, 
especially when combined with AO. 

The first observations of the GQ Lup companion have been 
carried out in K band in the night of Sept 16, 2005. The obser- 
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vations followed the standard scheme of sky nodding for back- 
ground subtraction. Eight nodding cycles with a integration time 
of 300s per frame where obtained. We have chosen the smallest 
pixel scale (12.5 x 25.0 mas) of the instrument, yielding a field 
of view of 0.8"x 0.8" to optimally sample the diffraction lim- 
ited core of the target PSF. The bright primary, GQ Lup A, was 
used as the AO guide star and was placed outside the FOV. The 
companion was centered in the field. 

The DIMV0 seeing during the science observations was 0.8- 
1.1". The strehl ratio, computed on a short integration of GQ 
Lup A was ~40% and more than 99% of the energy was encir- 
cled in a core of ~200 mas FWHM. 

H and J band observations have been obtained in the nights 
of April 24 and September 18, 2006, respectively. Ten (nine) 
target-sky nodding cycles of 300s exposure time have been taken 
in the H and J band, respectively. Five frames in each band had 
a sufficiently high signal-to-noise. 

The DIMM seeing for the five useful H band observations 
was 1.0-1.4". The strehl ratio, measured similarly as for the K 
band, was ~20% and more than 99% of the energy was encircled 
in a core of ~250 mas FWHM. 

No PSF calibration was possible for the J band since no ob- 
servation of GQ Lup A was taken. In situ strehl computations on 
the companion are not reliable because of the background con- 
tamination from GQ Lup A. We therefore judge from the seeing 
during the J band observations (0.8-1.1"), the airmass and the 
performance of MACAO on a strehl of 5-10%. The core size of 
the 99% quartile of the encircled energy is estimated to be of the 
order of 350 mas FWHM. 

2.2. Data reduction 

The data format of SINFONI is highly complex. The FOV of the 
instrument is sliced into 32 slitlets and remapped onto a pseudo- 
slit and dispersed by a grating. The original FOV has to be recon- 
structed from this spectrum. The success of this process depends 
on the proper identification of the slitlet positions on the detector 
and the correction of nonlinear distortions. 

The first step of the data reduction was thus performed by us- 
ing th e SINFON I data r eduction pipeline version 1.3 offered by 
ESO dJung et al.L 120 06). The reduction routines o f this pipeline 
where developed by the SINFONI consortium dAbuter et al.L 
2006). In this first step the images are cleaned from bad pix- 
els, flat fielded and wavelength calibrated. A distortion correc- 
tion is applied and the the FOV of SINFONI is reconstructed. 
The output of this procedure is a 3D fits cube containing about 
2000 images of the source in wavelength steps according to the 
chosen setting of the grating. 

A mean along the wavelength axis of such a cube is an equiv- 
alent of a broad-band image in J, H, or K band respectively. In 
Fig.[T]we show such images for all three bands. The companion 
is well visible in the centre of each image. Depending on the dif- 
ferent strehl ratios in each band, the target is still contaminated 
with light from the halo of GQ Lup A, whose core is outside 
the FOV. This contamination is spatially variable since the Airy 
pattern of the PSF of the primary is wavelength depended, even 
within a given band. This effect is strongest in the K band. 

For a proper extraction of the target spectrum we have to 
eliminate this contamination as far as possible. We dec ided to 
use the Starfinder package of IDL ( Diol aiti et all l2000h for an 
empirical PSF fitting of the companion in each of the ~2000 im- 
ages of the three observed bands. For this purpose a template 

1 Seeing value in the optical (A ~500nm), measured at zenith. 




Fig. 1. From left to right: SINFONI J, H and K band cube of 
the GQ Lup companion, integrated over the the wavelength in 
each cube. Top row. before subtraction of the halo of GQ Lup A, 
bottom row: after subtraction. North is up and east is left. The 
FOV is 0.8"x 0.8" 



PSF is created from the observation of a standard star right be- 
fore or after the target observation. The PSF of the standard star 
is thereby spatially supersampled from the ~2000 individual im- 
ages of the source. The Starfinder algorithm determines the flux 
of the companion in each wavelength bin by fitting the previ- 
ously created template PSF to the target. We found that one PSF 
template for each band is sufficient to ensure high correlation 
values (usually > 0.9) in the fitting process, despite the fact that 
the PSF shape is variable over the covered band, as explained 
before. In an iterative process the smooth background from the 
PSF halo of GQ Lup A is subtracted as well. The lower column 
in Fig. Q] shows the images of the companion after this subtrac- 
tion. The remaining speckles in the halo of GQ Lup A can not be 
filtered by this process. Since we fit a well sampled PSF to the 
target, the contribution from these localized remnants is however 
neglectable. 

Eventually we obtained a spectrum of the GQ Lup compan- 
ion from each nodding cycle in each of the three covered bands. 
These spectra are corrected for telluric absorption by the divi- 
sion through the spectrum of an early type telluric standard star. 
The following standard stars where chosen at the observatory 
to match the airmass of the science observations: HIP087140 
(B9V) for the J band, HIP082652 (B3III) for the H band, and 
HIP93193 (B9V) for the K band. These stars are intrinsically 
featureless, apart from weak helium and strong hydrogen lines. 
These lines have been manually fitted and removed. The result- 
ing spectrum was then multiplied by a blackbody spectrum to 
correct for the continuum slope of each standard star. The effec- 
tive temperature of the blackbody was retrieved from the liter- 
ature to r e yy=10500K for the stars of spectral type B9V and 
r e// =16500K for the B3III type star, respectively. Note, that 
we are in the Rayleigh- Jeans regime of the standard stars SED. 
Hence, the steepness of the continuum slope is not particularly 
sensitive to small errors in the effective temperature. After this 
procedure the resulting spectra of the companion are essentially 
free of telluric absorption lines and corrected for the throughput 
of the spectrograph. 

The individual spectra from each nodding cycle are com- 
bined by a weighted mean. The weights are derived from the 
correlation factor of the PSF fitting. The three combined J, H 
and K band spectra are displayed in Fig. [2] The Nyquist sam- 
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pled spectral resolution is R~2500 in the J-band and R~4000 in 
H and K band. Note that the signal-to-noise ratio of the spectra 
is ~ 100 in the J band but only ~30 for in the H and K band. 
These values have been retrieved from the standard deviation in 
each spectral bin of the five (J and H band) and eight (K band) 
individual spectra per band which have been reduced separately. 
The rather low signal-to-noise values origin from the spectral 
undersampling (1.5 pix per resolution element). Sub-pixel wave- 
length shifts between the individual nodding cycles induce artifi- 
cial noise. However, most of the small-scale features seen in the 
spectra are not noise but unresolved absorption lines. 

3. Spectral synthesis 

We compared our observed spectra with the theoretical template 
spect ra from the G AIA-Cond model v2 .0 dBrott & Hauschildtl 

2005) , updated from lAllard et al.l (1200 ll) , with improved molec- 
ular dissociation constants, more dust species with opacities, 
spherical symmetry, and a mixing length parameter 2.0 H p . The 
grid spans a range of 1000... 3000 K in effective temperature 
and -0.5. ..6.0 dex in surface gravity. The metallicity is fixed at 
[Fe/H]=0.0. The stepsize of the grid is 100K in effective temper- 
ature and 0.5 dex in surface gravity. Where appropriate we lin- 
early interpolated the spectral grid to refine the fits to the mea- 
sured spectra. At around 2600-2700K (the final effective tem- 
perature we derived for the companion to GQ Lup - see below), 
the effect of dust on the spectra is still small (it starts to affect 
the spectra massively around 2400K). Furthermore, dusty mod- 
els of this class were not available , thus we feel justified to use 
the cond models. The synthetic spectra are computed for a res- 
olution of 200000. We convolved the spectra with a Gaussian 
kernel to reduce the resolution to the one of the SINFONI spec- 
tra. Finally the synthetic spectra where re-binned to match the 
spectral sampling of the measured spectra. For the entire grid 
of synthetic spectra a ^-algorithm was applied to determine the 
best-fitting model spectrum, see Fig. [2] 

3.1. Spectral features 
J band 

A water vapor band absorption longwards of 1.30 /mi is very 
well fitted by the models and gives strong constraint on the ef- 
fective temperature. The strong metal lines (K I doublets at 1 . 169 
& 1.178 yum and at 1.244 & 1.256 jjm as well as Al I features at 
1.313 & 1.315 fim) are very sensitive on the effective tempera- 
ture and the surface gravity. However, their relative linestrengths 
can not be modeled accurately since the model grid has a fixed 
metallicity at [Fe/H]=0.0. Hence these features have been left 
out of the fitting process. The most prominent feature is the Pa 
f3 emission line at 1 .282 fim. This line shows an inverse P-Cygni 
profile in GQ Lup A. We identify the emission as being most 
likely caused by accretion in both objects. 

The continuum of the J band shortwards of 1 .25/mi looks ex- 
ceptionally bumpy and is only marginally fitted by the models. 
The steep decrease in flux from 1.15 to 1.2 fim is also seen in 
the spectra of the young field L dwarf, 2MAS S JO 14 1-4633, and 
most likely due to VO and FeH absorption ( Kirkpa trick et all 

2006) . However, in our spectrum it falls in a region that is most 
subject to remaining contamination from GQ Lup A. We thus 
decided to restrict the fitting process to wavelength longer than 
1 .2 //m. Moreover, the increasing slope from 1 .20 to 1 .23 fim ob- 
served in the GQ Lup companion was also observed in 2MASS 
J0141-4633 and other young, low-gravity objects like KPNO- 



GAIA cond 2.0 model: 1^2650^ logg=3.7 

1.6- 
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Fig. 2. From top to bottom: SINFONI J, H and K band spectra 
of the GQ Lup companion. The best fitting synthetic spectra are 
overplotted with a relative offset of +0.5. The Nyquist sampled 
spectral resolution is R~2500 in the J-band and R~4000 in H 
and K band. Note the strong Pa-/3 emission line in the J-band. 



Tau 4 and G 196-3 B dMcGovern et alll2004l) and is due to FeH 
bands. It is also not yet reproduced by our GAIA-cond model 
and has been consequently disgarded from the fit. 
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Fig. 3. Results from the fit of GAIA cond synthetic model spec- 
tra to the SINFONI J, H and K band spectra of the GQ Lup 
companion. Best value marked with a cross. Contours relate to 
1, 2 and 3<x errors. 



H band 

The continuum slope of the H band exhibits a strong triangu- 
lar shape that is usually identified with low-gravity objects. The 
nature of this feature is suppos edly the H2 collision in duced ab- 
sorption (CIA) as discussed by Kirkpatrick et al.l (120061 see ref- 
erences therein). The H band is essentially free of strong ab- 
soption lines from molecules and metals. The constraints on the 
effective temperature and surface gravity is the lowest among the 
three bands. 

Kband 

The K band spectrum shows an excellent agreement between 
model and measurement, both, in the slope of the pseudo- 
continuum and the depth of the most prominent spectral features. 
Again, water vapor absorption at both edges of the spectrum is 
a good indicator for the effective temperature. The I2 CO band 
heads longward of 2.29 fim are also very well fitted. We iden- 
tify a weak Na I doublet at 2.201 //m which is however slightly 
better fitted by a log g determined from J and H band. We note a 
non-detection of Br y line at 2.166 /mi. This line is in emission 
in the spectra of GQ Lup A where it shows a slightly asymmetric 
profile. 

3.2. Determination ofT eS and logg 

The results of the fitting process are displayed in Fig. [3] for the 
combined x 1 minimisation in a log g—T e g plane. The derived 
uncertainties (displayed in Fig. [3] as contours of 1, 2 and 3<x) 
stand for the systematic errors. We used a Monte-Carlo simu- 
lation technique to infer the statistical errors on the determina- 
tion of the model parameters. The uncertainties from the noise 
in the spectra where found to be much lower than the uncertain- 
ties from fitting errors. This demonstrates that the slope of the 
continuum is the major contraint as long as the spectral lines are 
not well resolved. 

While the best fit for the J band alone is = 2650 K and 
log g = 3 .70 dex, the best fit for the H and K band spectrum have 
slightly different effective temperatures and surface gravities, as 
can be seen in Fig. [2] Since the J band spectrum has the steepest 
X 1 curve, it is dominating the final solution, which we determine 
to be r eff = 2650 ± 100 K and log g = 3.70 + 0.5 dex. 



4. Results and discussion 

To determine the mass of the GQ Lup companion independently 
from evolutionary models, we need its surface gravity and ra- 
dius. The latter can be derived from a p reviously determ i ned K 
band magnitude K s = 13.1 ± 0.2 mag (Neuhauser et al., 2005) 
and a bolometric correction. Based on our new results for the 
effective temperature, we have to slightly re vise the bolomet- 
ric correction used in | Neuhauser et all (12005b . From table 6 of 
iGolimowski et alj d2004l) we determine BC K = 3.00 + 0.07 mag, 
interpolating over the objects with comparable effective temper- 
ature. 

With these values, we derive a luminosity of log(L/L ) = 
-2.25 + 0.24, slightly higher than in lNeuhauser et al.l (2005), but 
consistent within the error margin. From luminosity and temper- 
ature, we can then estimate the radius of the companion to be 
R = 3.50*Jq3 Rjup- Gravity and radius then determine the mass, 
for which we get ~ 25 Mj up , however, with large uncertainty. 
The mass can be as low as ~ 4 Mj up for lower gravity or even 
as high as ~ 155Mj up given the large uncertainty in the surface 
gravity. 

At face value, th e mass does not seem to be better contrained 
now, compared with iNeuhauser et all d2005l) . who gave 1 to 42 
Jupiter masses. However, we think that due to the low resolution, 
low signal-to-noise, and limited spectral coverage the surface 
gr avity and effective temp erature may have been underestimated 
in Neuhauser et al. (2005). Even though our mass determination 
is independent from evolutionary models, it still depends on the 
atmospheric models. Until the GAIA model is tested and con- 
firmed in the very low-mass and young age regime, our results 
are to be seen uncertain. However, we argue that the accuracy of 
a mass determined from spectral synthesis should be higher than 
from evolutionary models since the latter should use synthetic 
atmospheric models itself as an input parameter. 

Stassun et al. (2006) found an eclipsing double-lined 
spectroscopic binary, 2MASS J05352184-0546085 (hereafter 
2M0535), consisting of two brown dwarfs in the Orion star 
forming region. Their masses could be determined dynamically, 
hence independently from evolutionary models and model atmo- 
spheres. This makes 2M0535 most suitable for a direct compar- 
ison with the GQ Lup companion. 

Note, that the more massive object in 2M0535 is, however, 
cooler, which is not expected from standard models (e.g. Lyon 
or Tucson groups). Stassun et al. (2006) also note a systematic 
uncertainty of +300 K in the absolute temperatures, and can give 
only a precise temperature difference or fraction. 

We note that the companion to GQ Lup has roughly the same 
age but is fain t er and smaller than both components in 2M0535 
dStassun et al.L 120061) . see Table 1 . The effective temperature is 
comparable with the cooler component in this double. Judging 
from the temperature and radius of the nearly coeval objects, 
the GQ Lup companion must be lower in mass than each of the 
components in 2M0535, i.e. < 36 Mj up . 

Finally we note two recent results on GQ Lup: 

Marois et al. (2006) fit the RIJHKL-band spectral energy dis- 
tribution of the GQ Lup companion with a GAIA model to ob- 
tain physical parameters. They find a radius of 0.38 ± 0.05 Rq, 
and an effective temperature of 2335 + 1 00 K. Hence, the author s 
confirm temperature and luminosity of Neuhau ser et alj d2005l) . 
Their modelfits are reported to be gravity insensitiv and a log g 
of 3-4 is assumed. They use the evolutionary models from the 
Lyon group to obtain a mass of 10 to 20 Jupitermasses. 

McElwain et al. (2006) obtained an integral field J and H- 
band spectrum with OSIRIS at Keck, which is lower in dynamic 
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range, resolution, S/N, and spectral coverage than ours. They 'GQ Lup' onpage[T] 

confirm the luminosity, and temperature from Neuhauser e t al.l '2MASS J0141-4633' onpage[3| 

d2006l) and obtain a mass of 10 to 40 M Jup . 'KPNO-Tau 4' on pageH 

Both authors use the models of the Lyon group. We argue 'G 196-3 B' onpagelj 
that these models are not a pplicable for objects up to ~ 10 Myrs, 'GQ Lup b' on pagey 
following the discussion in Chabri eret al.l (120051) . 

Based on our and otherwise published results, one still can- 
not decide, whether the companion to GQ Lup is a massive 
planet or a low-mass brown dwarf. 
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Table 1. Comparing the GQ Lup companion to 2M0535 A and 
B 



Parameter 


GQ Lup b/B (a) 


2M0535 A 


2M0535 B 


Teff [K] 


2650 ± 100 


2650 ± 100 (b) 


2790 ± 105 (b) 


log (L bo i/L B ) 


-2.25 ± 0.24 


-1.70 ± 0.11 


-1.74 + 0.11 


Age [Myr] 


0-2 


0-3 


0-3 


Radius [R Q ] 


36 +0 15 
-Mi 


0.669 ± 0.034 


0.511 ±0.026 


Mass [M Q ] 


0.025!";™ (c) 


0.054 ± 0.005 


0.034 ± 0.003 



Remarks: (a) GQ Lup b, if planet or planet candidate; GQ Lup B, if 
brown dwarf or star. 



(b) ±300 K systematic uncertainty (Stassun et al. 2006). 

(c) Regarding the upper limit: our companion is fainter and smaller 
than both 2M0535 A and B, hence lower in mass than the components 
of 2M0535. 



